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PHASE TURBULENT FLOW OF HYDROGEN IN TUBES WITH 
TEMPERATURE RATIOS TO 23 
by Maynard F. Taylor 
Lewis Research Center 

SUMMARY 

The existing methods of correlating and predicting single-phase turbulent heat- 
transfer coefficients with variable properties give coefficients which are sometimes in 
poor agreement with measured values. Each prediction method can be used only over a 
limited range of pressure, temperature, ratio of axial position to diameter, and ratio of 
surface to fluid bulk temperature ratio. 

Local heat-transfer coefficients for hydrogen from 10 investigations using symmet- 
rically heated straight tubes were used to determine a single equation that will predict 
heat- transfer coefficients over a wide range of conditions including surface to fluid bulk 
temperature ratios to 23. 

The correlation equation predicts heat-transfer coefficients with acceptable accu- 
racy over a wide range of conditions 

Nu b = 0.023 Re“- 8 Pr“- 4 ^ 

where Nu^ is the bulk Nusselt number, Re^ is the bulk Reynolds number, Pr^ is the 
bulk Prandtl number, T g is surface temperature, and is bulk temperature. This 
equation does not predict heat-transfer coefficients with acceptable accuracy in the near- 
critical pressure and temperature region where the inlet temperature is between 45° R 
(25° K) and the transposed critical temperature (the temperature at which the specific 
heat of a fluid at constant pressure reaches a maximum), and where the inlet pressure 

O 

is between the critical pressure and 530 psia (3. 65 MN/m ). The correlation equation 
also accurately predicted variable property heat-transfer coefficients for helium gas and 
nitrogen gas. 




INTRODUCTION 


Hydrogen is attractive for use as a propellant for both advanced chemical and nu- 
clear rockets because of its low molecular weight, which gives a high specific impulse 
(i oc VT/m), and because of its excellent heat- transfer properties. These propulsion 
system concepts produce extreme conditions which are a severe test for existing methods 
of predicting heat-transfer coefficients. 

Most of the experimental heat-transfer data for hydrogen was taken in a tube at su- 
percritical pressures near the critical temperature. There are some data for hydrogen 
flowing through tubes with surface temperatures near those encountered in a solid core 
nuclear rocket. At present, no data exist for fluid temperatures above 2800° R 
(1556° K) for flow through tubes. 

In this investigation, all available single-phase heat-transfer data for hydrogen 
flowing turbulently through symmetrically heated straight circular tubes, along with all 
helium and nitrogen data with large variations in the physical properties, were studied. 
The range of conditions covered by the hydrogen data examined in this investigation is 
shown in table I. The data do not cover all possible combinations of the various condi- 
tions (e. g. , large T g /T b at large x/D). 

A new correlating equation is presented which correlates and predicts single-phase 
forced- convection heat-transfer coefficients for turbulent flow through tubes over a much 

TABLE I. - COMPLETE RANGE OF EXPERIMENTAL CONDITIONS 


COVERED BY INVESTIGATORS USING HYDROGEN 


Conditions 

Minimum 

Maximum 

Ratio of 
maximum to 
minimum 

Ratio of distance from 
entrance of test sec- 
tion to inside diameter 
of test section, x/D 

2.0 

252 

126 

Ratio of surface to bulk 
temperature, T g /T b 

1 . 1 

27. 6 

25 

Inlet temperature, T\, 
°R; °K 

45; 25 

573; 318 

13 

Absolute static pressure, 
2 

p, psia; MN/m 

18; 0.124 

2500; 17.2 

140 

Heat flux, Q/S, 

Btu/(sec)(in. 2 ); MW/m 2 

0.036; 0.059 

27.6; 45.6 

767 

Bulk Reynolds number, 
Re b 

7 500 

13 800 000 

1840 

Surface temperature, T g 
°R; °K 

53; 29.5 

5600; 3110 

106 


2 



greater range of conditions than was previously possible. This single equation can be 
used in place of the several correlating equations, each of which is recommended for use 
over a limited range of conditions. 


SYMBOLS 

constant from ref. 11 used in eq. (12) 

Cg exponent of T g /T b 

Cp specific heat of gas at constant pressure 

D inside diameter of test section 

G mass flow rate per unit cross-sectional area 

h local heat- transfer coefficient 

h av average heat-transfer coefficient 

h asymtotic heat-transfer coefficient 

^on specific impulse 

k thermal conductivity of gas 

L length of test section 

m molecular weight of gas 

Nu Nusselt number, hD/k 

Pr Prandtl number, c p/k 

Jr 

p absolute static pressure 

Q rate of heat transfer to gas 

Re b bulk Reynolds number, GD/ p b 

Re^ modified film Reynolds number, PfV^D/ 

Re g modified surface Reynolds number, p g V b D/p g 

Re b ^ modified Reynolds number, Pq 4 V b D/ P q 4 
S heat-transfer area 

T temperature 

T* transposed critical temperature (temperature at which specific heat of fluid at 

constant pressure reaches a maximum) 

V velocity 
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X 


parameter for specifying reference temperature T x = X(T g - T b ) + 
x distance from entrance of test section 

ju absolute viscosity of gas 

v kinematic viscosity of gas, p/p 

p density of gas 

Subscripts: 
av average 

b bulk (when applied to properties, indicates evaluation at bulk temperature T b ) 
c critical 

cal calculated 

ex experimental 

f film (when applied to properties indicates evaluation at film temperature T f ) 
i inlet 

s surface (when applied to properties indicates evaluation at surface temperature 
T s> 

0. 4 reference (when applied to properties indicates evaluation at 
T 0. 4 = T b + 4 < T s - T b> 

DISCUSSION OF EXISTING CORRELATIONS 

For many years, heat-transfer coefficients were predicted by evaluating the physi- 
cal properties and density of the fluid at a reference temperature T x = X(T g - T b ) + T b . 

Humble, Lowdermilk, and Desmon (ref. 1) correlated average heat-transfer coeffi- 
cients for heating and cooling of air flowing through a smooth tube with wall tempera- 
tures to 3050° R (1695° K) and surface to fluid bulk temperature ratios (T /TV.) from 
0. 46 to 3. 5 using the following equation in which the physical properties are evaluated at 
the film reference temperature X = 0. 5: 


Nu f 


= 0. 034 Re®' 8 Pr 


0.4 

f 



1 


(1) 


Taylor and Kirchgessner (ref. 2) extended the range of wall temperature to 5900° R 
(3275° K) and the range of T g /T b to 3. 9 and measured both local and average heat- 
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transfer coefficients for helium flowing through a tungsten tube. They found that the av- 
erage heat-transfer data for helium along with some air data from reference 1 could 
best be correlated by 


Nu ay f = 0. 021 Re 8 * 8 Prj?* 4 



( 2 ) 


while the local heat-transfer coefficients for helium were correlated by 

Nu f = 0. 021 Re 8 * 8 Pr 8 * 4 (3) 

Thompson and Geery (ref. 3) attained Tg/T^ as great as 16. 5 by lowering the inlet 
temperature of the hydrogen to 55° R (30. 6° K). At supercritical pressures, they re- 
ported two modes or regimes of heat transfer: A lower heat-flux regime A for which 
they recommend as the correlating equation 


Nu b 


0. 028 Re 8 * 8 Pr 8 * 4 



and high heat flux regime B with the correlating equation 


( 4 ) 


Nu 


b 


- 0.0217 Re 8 * 8 Pr 8 * 4 



( 5 ) 


Unfortunately, there was considerable overlap of the two regimes, which makes the 
choice of whether equation (4) or (5) should be used exceedingly difficult. 

McCarthy and Wolf (ref. 4) report heat-transfer data for both helium and hydrogen 
with Tg/T^ up to 11. 1 and inlet temperatures from 135° to 540° R (75° to 300° K). 
They reported that their local heat-transfer coefficients near the exit of their test sec- 
tions were best correlated by 


Nu. = 0. 025 ReP* 8 Pr 8 * 4 (_J 

' T b> 


-0. 55 


( 6 ) 


but for some reason the following equation was recommended: 


5 



( 7 ) 


Nu = 0. 025 Re?’ 8 Pr®’ 4 

fc> fe o 

An analysis by Deissler and Presler (ref. 5) showed that, in general, the best value 
for X in the reference temperature concept is 0. 4, and the available data could be rep- 
resented reasonably well by 


Nu 


0.4 


Re 


0.75 

0.4 


31 


( 8 ) 


Szetela (ref. 6) reported a few heat-transfer coefficients for hydrogen and showed 
that equation (3) was not suitable for correlating that data. 

Weiland (ref. 7) presented local heat-transfer data for helium and hydrogen flowing 
through a tube with an L/D of 252 and showed that equation (3) worked for x/D > 30 
and 1. 1 < Tg/T^ < 2. 5; however, for 2. 5 < T g /T^ < 4. 5, there was an effect of x/D, 
and equation (3) worked only for x/D > 160. 

Taylor (ref. 8) reported heat-transfer coefficients for helium and hydrogen in a 
study of the effect of dissociation of hydrogen at tube surface temperature to 5600° R 
(3110° K). It appeared that any dissociation of hydrogen at the tube wall had less effect 
on the heat-transfer characteristics than does Tg/T^. Equation (3) did not correlate 
the data with T g /T^ greater than about 3. 5. 

Dalle Donne and Bowditch (ref. 9) measured heat- transfer coefficients for helium 
and air and recommended using the following equation for 26 < x/D < 166: 


Nu b = 0. 024 Re®' 8 Pr®' 4 | — |exp - 

T b> 


0. 29 + 0. 0056 


( 1)1 


( 9 ) 


Perkins and Worsoe-Schmidt (ref. 10) reported heat-transfer coefficients for nitrogen 
with Tg/T^ to 7. 35 and recommended the following correlation based on bulk tempera- 
ture 


r 



and based on wall temperature 


6 



Nu = 0. 023 Re?" 8 Pr 8 * 4 
s s s 



( 11 ) 


0 7 0 7 

They recommend that the entry parameter 1 + (x/D)“ u * ‘(T /T^) ' be dropped from 
equation (10) for x/D greater than 40 and from equation (11) for x/D greater than 24. 
Simoneau and Hendricks (ref. 11) through the use of a simplified equation 


h = 



( 12 ) 


suggested that the variation of the transport and thermodynamic properties does not in- 
fluence the computation of the heat-transfer coefficient greatly. The constant is 
0. 048 for hydrogen, 0. 020 for helium, 0. 0042 for air, and 0. 00385 for carbon dioxide. 

Taylor (ref. 12) measured heat-transfer coefficients for precooled helium and hy- 
drogen flowing through a tungsten tube with wall temperatures to 5300° R (2945° K) and 
Tg/Tjj to 8. The use of equations (3), (6), (9), and (12) were studied by using the data 
of both reference 8 and 12. Simultaneously, Miller and Taylor (ref. 13) investigated the 
use of equations (3), (6), (8), and (9) to correlate the experimental data of Weiland 
(ref. 7) which is for an L/D of 250 (which is in the range of L/D of a nuclear rocket 
reactor). For large x/D, the exponent of T g /T^ in equation (9) was influenced too 
strongly by x/D. It was modified and resulted in the following equation: 


Nu b 


= 0. 021 Re 8 ’ 8 Pr 8 - 4 



0. 29 + 0. 0019 



(13) 


The data of Weiland (ref. 7) and Taylor (refs. 8 and 12) were best correlated by equa- 
tion (13). 

Hess and Kunz (ref. 14) made a study of heat transfer to supercritical hydrogen and 
recommended the equation 


Nu f = 0. 0208 Re 8 * 8 Pr 8 - 4 [ 1 + 0. 01457 — 


(14) 


Hendricks, Simoneau, and Friedman (ref. 15) reported the heat-transfer character- 
istics for cryogenic hydrogen at pressures from 1000 to 2500 psia (6. 9x10 to 
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17. 5 MN/m ) and recommended the use of equation (3) for x/D from 19 to 29, 27 to 46, 
and 30 to 69 depending on the inside diameter of the test section. 

Miller, Seader, and Trebes (refs. 16 and 17) measured hydrogen heat-transfer co- 
efficients with high heat flux, T g /T b to 27. 6, and supercritical pressures. They found 
that a modified form of equation (13) 


0. 8n„0. 4 | 


Nu q 4 = 0. 0204 Reg |Pr"- 1(1 + 0. 00983 


(15) 


correlated their data best for x/D from 30 to 60. This range of x/D also limited the 
T g /T b to 12. 6 or less. 

Hendricks, Graham, Hsu, and Friedman (ref. 18) reported heat-transfer coeffi- 
cients for hydrogen with inlet temperatures as low as 45° R (25° K) at supercritical pres- 
sures. They reported that in the supercritical pressure region the gross trend of suc- 
cess with which the data could be correlated with equation (3) increased directly with 
pressure. In general, however, all the supercritical data were used to test equations 
(3), (6), and (14), and no systematic attempt was made to study the effect of inlet tem- 
perature and inlet pressure. 

The heat-transfer coefficients for hydrogen with heat fluxes to 27. 6 Btu per second 

o 

per square inch (45. 1 MW/m ) are presented in reference 19 with the recommendation 
that equation (14) be used with the addition of a term C^, a function of the coolant tem- 
perature. The quantity varies nonlinearly from 2. 0 at 50° R (28° K) to 0. 85 at 
85° R (47° K): 


Nu f = 0. 0208 C T Re?* 8 Pr?* 4 (1 + 0. 01457 ) 

V V 


(16) 


A comprehensive study (ref. 20) of the recommended correlation equations (3) to 
(6), (9), (14), (15), and (16) was made with the conclusion that none of the equations 
showed any evidence of being greatly superior to the other correlation equations. Equa- 
tion (16) was selected for film-temperature-based correlations, and equation (5) was 
selected for bulk-temperature-based correlations. 

Petukhov, Kirillov, and Maidanik (ref. 21) measured heat- transfer coefficients for 
nitrogen gas with Tg/T^ to 6 and correlated the data with 


Nu b 


= 0. 0212 Re£- 8 Pr 8 ' 4 



(17) 


8 



where C 2 is a nonlinear function of x/D. In a discussion of reference 21, Taylor pre- 
sented the same type of correlation equation with different values for the exponent C 2 . 

He pointed out that the highest T g /T b data available should be used to determine the 
exponent, since small ratios are not greatly affected by the powers to which they are 
raised. 

The range of temperature, pressure, T g /T b , and x/D for the correlation equations 
recommended by references 3, 4, 15, 16, 19, and 20 are shown in figures 1(a), 1(b), 

1(c), 1(d), 1(e), and 1(f), respectively. Using the temperature- entropy diagram is a 
convenient method of showing the location of measured inlet pressure and temperature in 
relation to the saturation lines and critical pressure and temperature. 


METHOD OF CORRELATION 

The possibility of success in correlating heat-transfer data by including the ratio of 
surface to fluid bulk temperature ratio raised to a power which is a function of x/D was 
reported by references 9, 12, 13, and 21. Unfortunately, in all cases, the correlations 
were applicable only for large values of x/D. Obviously, a correlation which would be 
valid over the complete length of a tube (including the entrance region), as well as for 
wide variations in the physical properties, would be extremely useful. 

Because of the great interest in hydrogen for use as a nozzle coolant and propellant, 
this gas has been the subject of a number of experimental investigations covering the ex- 
tremely wide range of conditions shown in table n. This range of conditions affords not 
only an excellent test of any correlation but, in addition, provides a sufficiently wide 
range of parameter variation that their effects are clearly shown. 

The procedure of this report was to use the conventional type of equation presented 
in references 9, 12, 13, and 21. 


Nu b = 0. 023 Re b * 8 Pr^* 2 (18) 

(where C 2 is a function of x/D) and to determine an exponent C 2 which would best 
correlate the heat-transfer coefficients presented in the investigations shown in table n. 

The data of each investigator were analyzed individually. All the heat-transfer data 
points reported in the 11 investigations listed in table n were used in the present inves- 
tigation. Each investigator’s data were analyzed by a digital computer program. At 
each x/D, the exponent of the surface to fluid temperature ratio was computed so that 
the heat-transfer coefficient calculated by equation (18) was equal to the experimental 
heat-transfer coefficient. 
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Inlet fluid temperature, 


100 1 


Region covered by recom- 
mended equation 



Geery (ref. 3); equation (4) where 
TJTjj ranges from 1. 1 to 8. 7 and 
x/D ranges from 10 to 30; equa- 
tion (5) where T s /T b ranges from 
3. 5 to 16. 5 and x/d ranges from 
10 to 30. 



Wolf (ref. 4); equation (6) where 
TJT b ranges from 1.5 to 9. 2 and 
x/D ranges from 15 to 50. 



(c) Correlation from Hendricks, 
Simoneau, and Friedman (ref. 15); 
equation (3) where T s /T b ranges from 
1.4 to 9. land x/D ranges from 19 
to 29, 27 to 46, and 30 to 69. 


80 \— 


60 — 


40 — 


20l — 





(d) Correlation from Miller, Seader, 
and Trebes (refs. 16 and 17); equa- 
tion (15) where T s /T b ranges from 
1.6 to 12.6 and x/D ranges from 
30 to 50. 


(e) Correlation from reference 19; 
equation (16) where T $ /T b ranges 
from 6. 1 to 21 and x/D ranges from 
from 6 to 34. 


(f) Correlation from reference 20; equa- 
tion (16) where T s /T b ranges from 
1.8 to 21 and x/D ranges from 5. 7 
to 148. 


Figure 1. - Range of hydrogen inlet temperature and inlet pressure covered by correlation equations recommended by various 
investigators. 
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TABLE n. - TEST CONDITIONS FOR VARIOUS SOURCES OF HYDROGEN DATA 




Investigators 

Ref- 

erence 

Ratio of distance 
from entrance of 
test section to 
inside diameter 
of test section, 
x/D 

Ratio of surface 
to bulk 
temperature, 

V T b 

Inlet tempera- 
ture 

Pressure 

p 

psia MN/m 

Heat input 
Btu MW 

Mass flow rate 
lb/sec kg/sec 

Wall tempera- 
ture 

Inside diameter of 
test section 

°R 

°K 

°R 

°K 

in. 

cm 

(in. ^)(sec) 

m 

Thompson 

3 

10. 3 to 30. 9 

1. 1 to 16. 5 

55 to 


682 to 

4. 70 to 


0. 226 to 



77 to 

42. 8 to 

0. 192 

0.488 

and Geery 





56.6 

1392 

9.60 


13. 1 



1594 

885 



McCarthy and 

4 

5. 8 to 50. 2 

1. 5 to 11. 1 

135 to 

75 to 

32 to 






830 to 

461 to 

. 194 

.493 

Wolf 





311 

1354 

9.32 

14.8 

24.2 



2240 

1245 

.430 

1.093 

Weiland 

7 

2.7 to 252 

1. 1 to 4. 5 

222 to 

123 to 

250 to 

1. 72 to 


0. 392 to 

(a) 

(a) 


294 to 

. 188 

.478 





548 

304 


6. 89 


4.9 




1278 



Taylor 

8 

11,6 to 73. 1 

1. 3 to 5. 6 

557 to 


41 to 


0. 79 to 

1. 29 to 



738 to 


. 116 

.295 





573 

318 

79 


3.3 

5.4 



5630 

3130 



Taylor 

12 

11. 6 to 73. 1 

1. 4 to 8. 0 

263 to 

146 to 

37 to 

0.26 to 


1. 37 to 




194 to 

. 115 

.292 





304 

169 

68 


2.4 

3. 92 




2945 



Hendricks, 

15 

3. 4 to 78. 2 

1. 5 to 11. 0 

53 to 

29.4 to 


6. 89 to 





114 to 

63. 3 to 

.2110 

.536 

Simoneau, 





167 


17.2 


16.3 


0. 18 


666 

.438 

1.113 

and Friedman 
















Miller, 

16, 17 

4. 7 to 47. 4 

1. 6 to 27. 6 

51 to 

28.3 to 

458 to 

3. 16 to 

1.37 to 

2. 24 to 


0. 073 to 


59. 4 to 

.211 

.536 

Seader, 




69 

38.3 

2486 

17. 1 

24 

39.2 



1730 


.214 

.544 

and Trebes 
















Hendricks, 

18 

4.9 to 114 

1. 1 to 15. 6 

45 to 

25 to 

220 to 

1. 52 to 

0. 32 to 

0. 52 to 



45 to 

25 to 

. 188 

.478 

Graham, Hsu, 




115 

63.9 


5. 52 

3.3 

5.4 



1311 

728 


1.288 

and Friedman 
















Aerojet-General 

19 

6.7 to 33.9 

6. 1 to 21.4 

53 to 

29. 4 to 

696 to 

4. 80 to 

6. 4 to 




430 to 

239 to 

.147 

.373 

Corp. 




69 

38.3 

1371 

9.45 

27.6 

45. 1 



1681 

934 



Gladden 

(b) 

2. 0 to 127 

1. 4 to 6. 3 

65 to 

36. 1 to 

18 to 


(a) 

(a) 




139 to 

.335 

.851 

and Watt 




170 

94. 5 

36 







417 




^ot given. 

^Unpublished data, Lewis Research Center. 









1.4 
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1.2 

1.0 

.8 

.6 

.4 

.2 

0 

-2 

-.4 

-.6 

-.8 

- 1.0 


Source 


Exponent of sur- 
face to bulk tem- 
perature ratio 
(T s /T b ) in equation 


Ref. 9 

Refs. 12, 13 

Ref. 21 

Present 

rep 


(9) 

(13) 

(17) 

( 20 ) 


Ratio of distance 
from entrance * 
of test section to 
inside diameter of - 
test section, 
x/D 

26 to 166 
14 to 240 
5 to 100 
3. 4 to 252 


_L 


ejjort 


0 20 40 60 80 100 


1 

120 


f 

140 


I 

160 


180 200 


l 

220 


240 


Dimensionless distance from entrance, x/D 


Figure 2. - Variation of exponent of ratio of surface to fluid bulk tempera- 
tures with axial distance from entrance of test section. 


The computed exponents determined for all the experimental data were curve fit by 
using first-, second-, third-, and fourth- order polynomials in D/x„ Equal weight was 
given to each investigator's data. The fourth-order polynomial gave the best fit but was 
only slightly better than the more easily used and remembered first-order fit: 

C 2 = 0. 57 - ii®® (19) 

x 

D 


A plot of C 2 calculated from equation (19) is shown in figure 2. Exponents of Dalle 
Donne and Bowditch (ref. 9); Miller and Taylor (ref. 13); and Petukhov, Kirillov, and 
Maidanik (ref. 21) are also shown for comparison. 

Equations (18) and (19) are combined to give 


Nu^ = 0. 023 Re^ 


0. 8—0. 4/ T s 


Pr 



exp - | 0. 57 



( 20 ) 
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TABLE m. - PERCENT OF LOCAL HEAT- TRANSFER COEFFICIENT THAT DEVIATES 


I 

i 


LESS THAN ±25 PERCENT FROM EXPERIMENTAL HEAT-TRANSFER 


COEFFICIENTS FOR VARIOUS INVESTIGATORS 


t 

Investigator 

Reference 

Number of 

Percent of h cal that 

\ 




data points 

deviates less than 

[ 




used 

±25 percent from h ex 


Region 1, 

hydrogen 



Thompson and Geery 


3 

154 

68 


McCarthy and Wolf 


4 

249 

86 


Weiland 


7 

96 

99 


Taylor 

8 

and 12 

208 

97 


Hendricks, Simoneau, and Friedman 


15 

934 

88 


Miller, Seader, and Trebes 

16 

and 17 

560 

73 


Hendricks, Graham, Hsu, and Friedman 


18 

816 

89 


Aerojet- General Corp. 


19 

232 

98 


Gladden and Watt 


(a) 

425 

89 





3674 Total 

87 Average 


Region 2, 

hydrogen 



Miller, Seader, and Trebes 

16 

and 17 

60 

42 


Hendricks, Graham, Hsu, and Friedman 


18 

888 

39 





948 Total 

40 Average 


Helium gas 




McCarthy and Wolf 


4 

111 

96 


Weiland 


7 

96 

99 


Taylor 

8 and 12 

152 

98 

< 




359 Total 

98 Average 


Nitrogen gas 




Perkins and Worsoe-Schmidt 


10 

88 

97 


Unpublished data, Lewis Research Center. 
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A computer program by Goldberg and Haferd (ref. 22) was used to calculate physical 
properties of hydrogen for use in equation (20). 


DISCUSSION OF RESULTS 


The measured heat-transfer coefficients of the investigations shown in table n were 
compared with the heat-transfer coefficients calculated by equation (20) for the experi- 
mental conditions. The percent deviation was calculated by 



100 = Percent deviation 


( 21 ) 


The total number of data points used from each investigation and the percentage of them 
which had deviations within ±25 percent are recorded in table in. The hydrogen data are 
separated into regions 1 and 2, as shown in figure 3, so that the effect of inlet tempera- 
ture, inlet pressure, and transposed critical temperature could be studied. The results 
shown in table m indicate that equation (20) can be used with confidence over the wide 
range of conditions covered by region 1. In region 1, 87 percent of the 3674 calculated 
heat-transfer coefficients deviated less than ±25 percent from the measured values. 
Region 2 defined by 45° R (2 5° K) < T^ < T* and P c < Pj < 530 psia (3. 65 MN/m^) is 
often referred to as the near- critical region. In region 2, only 40 percent of 948 calcu- 
lated heat-transfer coefficients deviated less than ±2 5 percent from the measured values. 
Equation (20) should be used with caution in this region. At present, however, there is 
no other equation which predicts heat-transfer coefficients with much more success. The 
inability of equation (20) to predict heat-transfer coefficients in this region may be a re- 
sult of the lack of accurate transport properties or to a mode of heat transfer which is 
not approximated by single-phase heat transfer. 

In addition to the hydrogen data, 88 nitrogen gas data points reported by Perkins and 
Worsoe-Schmidt (ref. 10) and 359 helium gas data points reported by McCarthy and Wolf 
(ref. 4), Weiland (ref. 7), and Taylor (refs. 8 and 12) were predicted by using equa- 
tion (20) with 98 percent of the calculated heat-transfer coefficients deviating less than 
±25 percent from the measured heat-transfer coefficients. 

A fair comparison between the correlation equations (recommended by the investi- 
gators listed in table n) and equation (20) could not be made since these equations were 
not recommended for use at small values of x/D, where the largest Tg/T^ usually 
occurs. A comparison could be made by using both the data from reference 19 with 
x/D from 6. 7 to 33. 9 and the recommended correlation, equation (16), from refer- 
ence 20 for x/D from 5. 7 to 148. Figure 4 shows the results of comparing equa- 
tion (16) of reference 20 with equation (20) of the present investigation. As shown in fig- 
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(a) Correlation using equation (16) recommended by (b) Correlation using equation (20) recom- 

references 19 and 20. mended by present investigation. 

Figure 4. - Variation of deviation of the calculated heat-transfer coefficient from measured heat-transfer coefficient with 
local fluid bulk temperature (232 data points from ref. 19). Ratio of distance from entrance of test section to inside diameter 
of test section, x/D, 6. 7 to 33. 9. 

ure 4(a), 66 percent of the heat-transfer coefficients predicted by equation (16) deviated 
less than ±20 percent from measured values, and 40 percent deviated less than 10 per- 
cent. Figure 4(b) shows the results of using equation (20) to predict heat-transfer coef- 
ficients: 94 percent of the calculated heat-transfer coefficients deviated less than ±20 
percent, and 64 percent deviated less than ±10 percent from the measured heat-transfer 
coefficients. 

Equation (20) has predicted heat-transfer coefficients for hydrogen which are in 
good agreement with experimental values over a wide range of conditions including 
values of x/D from 2. 0 to 252. However, the wide variation of the exponent Cg with 
x/D for values of x/D less than about 5 indicates a need for caution in using equa- 
tion (20) in this range. None of the investigations listed in table II were made to study 
heat transfer in the entrance region. Those with large Tg/T^ usually have large sur- 
face temperature gradients which are not conducive to accurate heat-transfer data in the 
entrance region. 


I 
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Figure 5. - Variation of local to asymptotic heat-transfer coefficient for 
water with ratio of distance from entrance of test section to inside diameter of 
test section. 


The ratio of the local heat-transfer coefficient h to the asymptotic heat-transfer 
coefficient h g for water for several values of x/D from 0„ 19 to 6. 3 and a Tg/T^ of 
1. 06 were taken from an experimental investigation by Harnett (ref. 23) and are shown 
in figure 5. The analytical curves of Deissler (ref. 24) and Sparrow, Hallman, and 
Siegel (ref. 2 5) for a Prandtl number of 10 and constant properties are shown with the 
curve predicted by equation (20) for a Tg/T^ of 1. 06. For this particular case with 
T g /T b near 1, there is good agreement between the experimental data, the analytical 
results, and the prediction of equation (20). Equation (20) was also used to predict heat- 
transfer coefficients for air with T g /T^ from 1. 10 to 1. 91, as presented by Magee 
(ref. 26). In figure 6, the ratio of the measured local heat transfer to the asymptotic 
heat-transfer coefficient is compared with curves predicted by equation (20) using the 
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Figure 6. - Variation of local to fully developed heat -transfer coefficient for air with ratio of distance from 
/ entrance of test section to inside diameter of test section. Prandt! number, 0. 71. 

local experimental values for T g /T b , x/D, and Reynolds number showing an effect of 
both T g /T b and Reynolds number. 

The need for caution in using heat-transfer equations for predicting wall tempera- 
tures has been demonstrated by Miller and Taylor (ref. 13). Since equation (20) is em- 
pirical, it is recommended that it not be extrapolated beyond the conditions for which it 
has been tested. 


CONCLUSIONS 

Variable-property single-phase hydrogen heat-transfer data from 10 investigations 
that used symmetrically heated straight tubes have been correlated over a wide range of 
conditions by the single equation 



where Nu b is the bulk Nusselt number, Re b is the bulk Reynolds number, Pr b is the 
bulk Prandtl number, T g is surface temperature, and T b is bulk temperature. The 
same equation was also used to correlate helium and nitrogen data. 

The range of conditions over which the correlation applies for hydrogen is 
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1. 1 to 23 


Ratio of surface to bulk temperature, Tg/T^ . „ , 

Ratio of distance from entrance of test section to 

inside diameter of test section, x/D 2. 0 to 252 

Heat flux Q/S, Btu/(sec)(in. 2 ); MW/m 2 0. 036 to 27. 6; 0.059 to 45. 7 

Bulk Reynolds number, Re^ 7 500 to 13 800 000 

Surface temperature, T , °R; °K 114 to 5630; 63 to 3130 

All possible combinations of these conditions are not available. 

The hydrogen heat-transfer conditions were divided into regions 1 and 2 in the inlet 
temperature and inlet pressure plane. The deviation of predicted heat-transfer coeffi- 
cients from measured values in each region is as follows: 

1. In region 1, 87 percent of the 3674 calculated heat-transfer coefficients deviated 
less than ±25 percent from the measured heat-transfer coefficients. The correlation 
equation is recommended to predict heat-transfer coefficients in this region. 

2. In region 2, 40 percent of the 948 calculated heat-transfer coefficients deviated 
less than ±2 5 percent from the measured heat- transfer coefficients. Since no other 
equation can predict heat-transfer coefficients much more successfully, the correlation 
equation is suggested for use with caution in this region. 

The correlation equation was used to predict heat-transfer coefficients for helium 
gas and nitrogen gas. Of 447 heat-transfer coefficients that were calculated, 98 percent 
deviated less than ±2 5 percent from the measured heat-transfer coefficients. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 27, 1967, 

122-28-02-16-22. 
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